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THE MAGNETIC STATE OF THE EARTH AT EPOCH 1885.0 


Adolph Schmidt' 

OoHm 


The general theoretical developments on which the investiga- 
tions presented here are based, were presented in an earlier 
issue of this journal (XII, 2, 1889). A report of the most 
important results proviaed by the application of these develop- 
ments to the state of the earth's magnetism in the year 1885 has 
also been published ( Abhondlungen der k. bayer. Akodemfe der 
Whscnsch., II.Kl., XIX Vol . J Abth) On the following pages the 
fundamentals and the results of this investigation will be presented 
in greater detail and in addition, the calculations themselves 
will be presented. This is done with the intent of creating a 
reliable and comfortable basis for future investigations of a 
similar nature, but not because the results themselves might have 
a conclusive significance. These results doubtless require 
significant improvement and it has long been my intention to 
perform a definitive recalculation for this period of time as soon 
as the needed information becomes available. The presentation 
here is based entirely on the values of the earth's magnetic force 
components at 1800 points of the earth's surface derived by Dr. 
Neumayer. The observations on which these values are based, extend 
back to about 1887; the vast majority come from the time before 
1885, to which the derived values presented in the earth atlas 
of magnetism pertain. The determination of these latter values 
had to be by extrapolation in most cases; this necessarily reduces 
their validity. This circumstance was unavoidable since the atlas 
naturally was to provide a representation for a short time segment. 
But in addition to this problem, there is the deficiency of the 
observation material for wide regions, as Dr. Neumayer discussed 
in detail in the notes on his atlas. It is clear that progress 
has been achieved in two ways through the incorporation of recent 
observations — which it is hoped will be expanded considerably in 
coming years. The scope of valuable material has increased, both 
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in uniformity of geographic distribution, and the use of 
observations symmetrical to the normal epoch increases the 
reliability of application to the earth. If these considerations 
should make a future repitition of the present calculations seem 
expedient, then in addition it should be noted that observations 
made above 60° N-latitute have been entirely excluded. 

In spite of the described, generally unavoidable deficiencies, 
it is hoped that the reporting of provisional results will not be 
thought unjustified, not only because of the simplification this 
means for a final working, but also because the 
anticipated observations from the South Polar regions will 
have to be delayed for several years, and also because no 
signficantly better results are likely to be available for some 
time . 


The two papers mentioned above — which will be referenced 
below as A and B, contain such a detailed presentation of every- 
thing not relating exclusively to the performance of the numerical 
calculation, that I can limit this discussion almost entirely 
to an exposition of these calculations. Thus, rep .titions have 
been prevented, except where absolutely necessary for the 
cohesion of this presentation. 

Survey of Mathematic Aids in the Expansion 

The empirical basis of the entire investigation is formed by 
the maps of the geomagnetic elements E,6. i constructed by Dr. 

Neumayer for the beginning of the year 1885, or rather, by the 
values of these quantities taken by him for 1800 points where the 
meridians of 0°, 5°, 10°. ...355° East longitude from Greenwich and 
the parallel circles of 0°, 5°, 10°.... 60° North and Sough geographic 
latitude intersect. The mentioned, detailed text of the atlas of 
geomagnetism (the 4th part of Berghaus' Physical Atlas) provides 
information about the materials used in construction of the maps 
and about the applied methods of map-making; it is thus unnecessary 
to discuss these matters any further here. 
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From Che values of the elements H, i, i , those of the com- 
ponents X, Y, Z were derived. These latter are presented in 
table III. ’X’ means North, 'Y' means East, '2' means the 
downward, positive-measured component of force, so that the 
arrangement of positive semi-axes agrees with the present standard. 
The unit of measure used here and in all numbers in the present 
report, is 0.1^ cm"~^g^s~^, i.e. the unit of the last place 
which still has some relevance in variation observations. 

Prof. Eschenhagen suggested the designation y as a remembrance 
of Gauss (see Terrestrial Magnetism, Vol. I, p. 57, note 2). I 
will use this designation hereafter. 

The values of the components were then presented on each 
latitude by means of trigonometric series as functions of geographic 
longitude 4. The coefficients of this series developed to 4th 
order terms are presented in table IVa, b, c. They formed the 
starting data for my own calculation. 


X = Iti ro» i+X| «n2+ + Xi«m4i 

r 7, coff i+ ii tin 4* WB 4i 

Z — rof z+ifi «n i+ +i4««4i 


If the flattening of the earth is ignored, then the obtained 
numbers could be expressed by spherical functions of geographic 
latitude. But this is directly possible only for Z, since X and Y 
are inconstant at the poles because they approach equivocal 
expressions of the form 

ncota^a,) c,ein(i—ai) at the North pole 

—rtcota—nj) and at the South pole. 

So whereas Z can be developed without change, X and Y must be 
represented by expressions formed in connection with Z, or by 
themselves, and these expressions must be free of all discontinuity. 
This can be done in a variety of ways. A limitation is introduced 
by the requirement that the selected expressions should permit a 
simple and a closed series development leading to a derivation of 
the potential on the earth's surface. The simplest possible 
values in this case, which are sufficient for a unique definition 
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of the force vector, are X sin u, Y sin u and Z, if u denotes 
the complement of the geographic latitude. Besides these, the 
following might also be taken into consideration: 

rcottintti, XMinX—Tatufctl, Z 

and —Xmueotl-^Ttinl^Ztinut^l, —XcMutinX+Teetk—ZnnumnX, Z»inH~~Zec$H. 

The second group represents the components of force in three 
fixed axes, i.e. rectified at all points; these axes are parallel 
to the earth radii to the equatorial points of 0° and 90° East 
longitude and to the North pole. 

Now if the deviation of the earth's surface from the 

spherical is to be taken into account, as is the case here, then 

this necessitates a modification of the calculation (see A, p. 13 

B, p. 4). First, the geographic latitude has to be replaced by 

the geocentric latitude. Its complement, called », is defined 

by the equation: .. 

tffp =a VT+***9*t = [0.0014542] 

where the bracketed figure is the usual abbreviation foi num log. 

The value of ** used here is 0.00671922; it corresponds to the 

Bessel factor for flattening, 1:299.1528. The computed values 

of p belonging to u = 0°, 5°, 10°...90° have been rounded off to 

whole seconds of degrees (see B, p. 5): 

0*0' 0* 5*1' 0* 10*l'a«' 16*2'58* 20*8' 42* 25*4"25' a0*4'58* 85*6'25' 40*6' 40' 

45*6'46* 80*5'40' 55'8'24' 60*4'89' 65*4'24' 70*8'42' 75*2'82* 80*1'68’ 85*1' 0* 

90*0' 0*. 


For Uj^ = 180° - u, we have = 180 “— p. All other calculations 
are based on the rounded values given here, and not on the equa- 
tion presented above. 

Another deviation of the calculation from that of a sphere 
is that instead of the force components X, Y, Z, we have to use 
the slightly different quantities «x, 4F, where: 
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a = yi*f «*ea»r* 

tinp 
* tm u 


A^ylT? 


eotc 

* CWII 


The quantities to be developed according to spherical functions 
of the argument p are now: 

aXmo 0Ttmp fZ. 


For the parallel circle of u = 0®, 5°, 10°. . . .90°, , . .180°, 
the logarithms of the coefficients contained herein are (their 
numerical values have been given in B, p. A7, table II, in addi- 
tion to those for «, 4 and r ) : 


tt a iinv 

logfirinp 

lofr 

1 

i 

ti iloffu$mp 

%4«nv 

lagr 


0* i — oo 

— oo 

0.0000000 

180* 

45* ' 9.8509361 

9.8516644 

9.9992717 

185* 

& : 8J4S1807 

8.9431918 

9.9999889 

175 

50 1 9.8854533 

9.88^78 

9.9991455 

ISO 

10 '9.34S4871 

9.9425311 

9.9999560 

170 

55 9.9143187 

9.9152956 

9.9990231 

125 

15 '9.4157098 

9.4158075 

9.9999023 

165 

60 ; 9.9382561 

9.9393477 

9.9989084 

120 

SO 9.5366174 

9.5367880 

9.9998294 

160 

65 9.9577935 

9.9589887 

9.9988048 

115 

95 19.6283366 

9.6285970 

9.9997396 

155 

70 9.9733253 

9.9746099 

9.9987154 

no 

30 ; 9.7011483 

9.7015128 

9.9996355 

150 

75 f 9.9851378 

9.9864949 

9.9986429 

105 

35 j 9.7605416 

9.7610211 

9.9995205 

145 

80 P 9.9934389 

9.9948494 

9.9985895 

100 

40 ' 9.8097714 

9.8103734 

9.9993980 

140 

85 1 9.9988663 

9.9998095 

9.9985568 

95 

45 [ 9.8509361 

9.8516644 

9.9992717 

135 

90 1 0.0000000 

0.0014542 

9.9985458 

» 


Since X, Y and Z have already been developed by i, the only 
problem left to solve is the representation of the coefficients 

ak^tinp, aKmtinp\ fiLmtinp-, 

by spherical functions of m-th rank (PZ, PV' ) • For each of 

these coefficients, 25 values are known which belong to the 
parallel circles of geographic North-pole distances u = 30°, 35°.. 
..150°, which I will denote by a second, lower index i = 1,2,... 25. 
For i+rnM, we have *v+rr *= 180*, thus sin v^ = sin v^, and etc. 

Now it is clear that the sums: 

mkrn,i>inVi’i‘afkm,fPMPc, JH^r 

£ g ^ 

are even functions of cos », and that the corresponding differences: 



are uneven functions of cos v . The former depend only on spherical 

functions pI! with even difference (n-m) of the two indices; the 
m 
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latter depend on those with uneven difference (n-m), and since 
the observation data can be expressed completely by those sums 
and differences (due to their symmetrical distribution about the 
equator), the unknowns — the coefficients of spherical functions — 
break down into two separately determined groups. With regard to 
this circumstance — which considerably simplifies the numerical 
calculation — I do not intend to report the quantities . . .. 

(in table Va, b, c), rather only the cited sums and differences. 


The spherical function (or in Gaussian notation) 

is defined by the equation: 

. (f>— m) (w— fw— 1) (n— ffl— 2) rw—m— 8) , _ . 

'*■ 2.4.(2fi-l)(2»— 8) 00# 

For the functions up to 7th order used below, i.e. for those in 

which n^7 , I intend to compile the resultant series since it 

is convenient for many purposes to have the numerical values of 

the coefficients at hand. As abbreviation, I write instead 

m 

of , furthermore, c instead of eotv and s instead of $inv . 



PJ«=c g. 

p; = . P* « #e 







Pi- 


i1- 


•’('•-T') 


PJ 

PJ-f* 


£ 
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PJ « c»-ge*+|gc»-^ 





Pi - 

'i — (‘•-n"+ns) 

n - .•(.•-i) 



1 

1 

p;-rf* 

p; ««•« 

p; - 


For a numerical calculation of the functional values and 
in particular of their most-frequently used logarithms, a repre- 
sentation by products is preferred to this one using sums. 

p 2 is a whole function of mv* , except for a factor mt* 
m 

or If we break this down into its real, linear factors, 

we then obtain: 

JPZt(eoiv) ■* «i)(m«v— « i) («w«— 

I ••••• 


In this form the equation is valid for all cases. Now if 

n-m is uneven, then the independent value ■■ o appears and 

obtains the factor cc$p. , as it must be. 
m 


Now if we set: «, ■■ 

and note that: (eo$ v ^ eo$ a) ieo$v-i- cm a) — — «m(v+«)nN(v-«) » 

then we find: 


with 

and for uneven n-m: 


> ^ w * 180^, * 180^..... 

“ 90^. 


In order to make this convenient formula applicable for the 
numerical calculation, it is sufficient to compute the occurring, 
constant angles «t, «ii • . «n-«i for the various spherical functions 
one time only. I have done this and present the results in the 
overview below. Naturally it was expedient to carry the computa- 
tion out so that the results will be useful for all future appll- 
cations. Therefore, although a much less stringent computation 
would be sufficient for present purposes, the calculation was 
performed with 10-place logarithms (from Vega-s Thesaurus) and 
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the results are presented to 4 decimal places of seconds of arc. 
Their uncertainty is approximately 0”.0002. It hardly need be 
stated that in the representation of P”, the roots of the equa*> 
tion Pq = 0 already computed by Gauss could be used. The appro- 
priate angles ut are found in the tables published by Prof. Seeliger 
for the Neumann method of coefficient computation of spherical 
function series; the results are accurate to tenths of an arc- 

second ( Sfhcungiberichte der meth.-phyi. d. h. Aioidemtt d. Wtisenichoft 
zu Muenchen, 1890. poge499.) 


«, — M*44' 8T1971 
Oi ■» l25M5'6ir80» 

p, at mm 63*86' 5rS156 
aj • n6*88'64nW2 

p, a, 49* 6'28rT792 
-*»• 0, — 180* 81' 86:2906 

a, a* 28* 1' 2T4282 
a, — 154*6r57T8768 

ai 40* 8'iri09l 

* ‘ 04 wm 189* 84' 42:8909 

a, — 84*44' S:i97l 
•' ai I28*18'8i:8029 

- a, — 21* 10' 86:8448 
04 mm 188*49'28:i888 

ai i« 88*82'4i:?201 

* ' ai f 146* 7' 18:2799 

0, K 48*89*84:7020 
04 184* 0*28:2980 

pt ot mt, 86*81' 4:2482 
■'•‘‘a, * 121*28*88:7648 
a, - 16*21*28:2940 
181*88'8i:7060 

„ a, V 29*20*18:6861 
ai a 180*89*41:8689 
, a, . 89*41*41:9908 
atom 140*18*16:0098 

, a, M 60* 9*8e:9l08 
04 i> 129*80*28:0892 

^ a, . 61*17*22:1467 
118*42*87:8888 


JC («• r) a* MR fftin (t>4*a.) 


P ». 

•• 


a, — 89*18'8r478? 
a» >a 140*46* 6:8^ 

a, IK 80*88*20:i802 

• 04 mm 149*26*89:8698 


a|i« 90* 

0i « 70* r27:4iii 
ai >• 109*82*82:68^ 


ai » 90* 


PJ: •* " 

* as — 


6r4r8r4446 
112* ir 27:8884 


as aa 87*28*18:8042 
04 » 122*84*46:1968 
as — 78*26*88:8884 

o, ~ 106*84*81:6766 


-a 80* 


as — 90* 


jp». *1 “ 

^ al- 


as « 48*86*2ri779 
Ml — 181*88*81:8221 

as «> 62* 2*28:4878 
ai 1- nr8rM:84t8 

a, V 78*89*81:0887 
as » 104*80*88:9478 


70*81*4r6067 

109*tri6:89tt 

7rii* 41:7914 
108*4ri8»M 


• 90* 


90* 


JH . •• « 
* «1 ■■ 


a,-. 42* 8*16:7888 
o, mo 187*81*48:2467 

as -a 88*48*20:0984 
as — 126* 16' 89:9046 

as « 68* 6 * 8 r 6168 

04 mm 114*88*82:4888 

a, « 76*88*89:8898 

a. » 108* 4* 0:6402 


90* 


^4- 


at 

•1 


72*2T 8:7878 
107*88*64:8428 

66* 8' 21 : 341 6 
118*86*88:7884 

7r88* rtm 

108* 4*82:6861 

9cr 


78*88*88:8908 
106* 6* 7:6098 


M eo* 
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In this overview the functions p” and P” - have been omitted 

n n-i 

because they already appear in the original formulas as products 
and , 


The various spherical functions Pjj do not differ considerably 
from each other's average values. This is a disadvantage for 
numerical expansions which is noticed all the more, the farther 
the series is carried. Therefore, I have added to the functions 
p” (see B, p. 6, 7) constant factors r” of such magnitude that 
the quadratic average of the product rl-p”, which I call taken 
over the entire spherical surface, is equal to 1 for all values of 
m and n. For this purpose we must set: 


1.8.8.. . (8s with ^ * 1. n « .... »t 

(this results from the known properties of the functions P*^). 

m 

For the factors of the 7th order functions using this formula, we 
obtain the following expressions, whose values have been reported 


In B, p. 47, 

table III 

to 

7 significant 

figures . 




tt 

1 

rj 

- jKT 



108 
rJ - -g- 

r* - f HI 

»• — Ig 

^5- 

^Hf 
» ' 

r{ - /T 

rj 


•1- 


r* wm 

r» jllis 





fj 

- jll! 

rJ- 

|doi 


1 

i 

fS-|K5TI5 






•1- 

•yy55 

fj . Ay?o 

rJ « 

rS •> jjiCTi 




— vir r* 


r* — f* — Y*®® *** 

ft - ft m, rl 

I# * K 


M 


r} 


The logarithms of the functions Jlji} . ^ for the values of • 
coming into consideration here, are found in table I; the function 
values themselves are presented in B, p. 48/50, table IV. The 
reported numbers (in whose computation I have not yet applied 
the products stated for P|JJ) have been calculated with 7-place 


9 


ORIG^L PAGE m 
OF POOR QUALITY 


The logarithms of these numbers are:* 


M;e:' 

0 

1 

t 

1 • 

4 

s 

• 

7 

0 

0.0000000 

0.SSIM07 

0.ftt557«S 

0.St04t»0 

Uliotst 

141M4II 

2.71444S7 

44144440 

1 


0.SMS407 

0.m04ft7 

1 0.90MS47 

l.SSOISfl 

i.4f7ftrii 

IJStOOM 

11474040 

t 



0.II70IA7 

0.70fM47 

1.0«f444< 

Mossiest 

1.7I144M 

44444444 

s 




04S04Rf0 

0.7mm 

USStITt 

' 14444440 

14444414 

4 





0J4604M 

04467417 1 

14UI144 

14744740 

• 






04447417 ' 

1 04147444 

14744440 

• 







04441447 

0.4741444 

f 






i 

1 


04H1444 


logs while omitting the last digit In the key values attainable 
only through the use of multi-place tables. This place will thus 
sometimes be Inaccurate by somewhat more than half a unit; In the 
7-place values, sometimes even by one whole un5.t. 


The table of functions r" Is Indeed sufficient for the 

m 

derivation of rJf«i for the values of r contained therein; 

but It is convenient In many regards to be able to determine the 

coefficients k . . . directly from the applicable series for 
m ni ^ 

0Trinv,r2, In order to do this, tables of the logs of 

are needed. Since these tables can be used 
repeatedly, I have prepared them for Inclusion (as table II) at 
the end of the work. It seemed sufficient to cite values of these 
figures rounded to 4 decimal places since In future calculations 
of potential, only the deviations from the values determined here 
will come Into consideration — that Is, relatively small values. 

In this severe rounding, and loff(Sl:0nnt) differ at most 

by 15 units In the last place. Therefore, X have reported only 
the compilation of values of the former function and specified 
the (dependent on the angle v. Identical with ) difference 

(Iv unfitly tin 0 ) , which Is to be subtracted from It In order to 
obtain The pile-up of rounding errors occurring In 

many numbers, which could be eliminated by addition of -fl or -1 
to the last decimal, Is of no Importance to the purpose of the 
table. I thus felt justified In omitting any reference to this 

*Llne 3, 4 ami 5. In log r^, log r^, log rf & log r| , In the last 
place write a "6” instead of a ”7”. 


ORIGINAL RAGE 13 
OF POOR QUALITY 

problem In Che cable 

Those figures under Che funcClons (M:«Wnr) end (Ri,:0tmt) 
whose lower index m Is equal Co zerci beccxne InflnlCe for r o 
and l.e. aC boCh poles. Thus, Che pare of Che expansion 

of X and Y dependenc on Chem musC be .:ransformed, and ic is 
expedlenC Co use Chis Crans formed on for Che ocher values of r. 

The expansion for X follows easily from Che ocher reporCed 
condiCional equaCions for Che coefficiencs of Che spherical funeCion 
series, according Co Che infiniCe expressions: 

Ji'-ywTJ! and J8r‘-r^4.+l)lS 

• fm« mtmw 

where only m is Co be replaced by # for Chis represenCaCion (see 
also B, p. 25, 26). Thus, Che logs of Chese values have been 
incorporaced inCo Cable Ita. 

The coefflcienCs of Che series used in Che represenCaCion of 
mitinv, 0Ttmv and are subjecC Co cerCt’.n conditions, some of which 
have already been discussed. Since Chese wlH be discussed in 
deCail in the coming sections, this passing mention will suffice 
at this time. 

The expansion of the force components in fixed directions is 
simpler in an analytical respect — since only one very simple 
conditional equation has to be taken into account — it was already 
discussed above (p. 3). Therefore, it may be permissible to go 
into this in brief, even though no use is Co be made of Che 
appertlnent series. I will call Chose components x, M, Z and 
write Chem in the following form: 

X me Xmmu—M.mu 
Xtmu w 

X w —iX.$mum l-k-B mu 


from which follows: 
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These formulas are generally valid, but can only be used for 
a sphere since the expansion for the ellipsoid is performed by 
r and not by u. From this we see that for s, h, z we obtain 
limited spherical function series when values for Xiinu, Trinu^ Z 
are given, and that under certain conditions — the ones mentioned 
iust above — the reverse will apply. In order to perform the 
real transformation, a number of identities has to be used through 
which the products 

Si^(cMu).dgu 

are converted into sums of spherical functions, to which, however, 
for the last two products, expressions of another form appear, 
which consequently enable the conditional equations in s and H 
applicable for X sin u and Y sin u to disappear. As can be 
derived from the known properties of spherical functions, we have 
an abbreviated notation with J, — 


Tm 


rill *— I ri 

__ **111 nB4-l tm *'■»♦> 


for iii>0 


These formulas generally apply if we specify that r^ 
is zero for n<w. From this follows the recursive formulas below, 
which provide the solution for the last two problems: 


ja . «»«•« = ja-». 

r;.J 2 (ri.1)* 

rl 


•«« ’Comu 


ri>, H + m rl 


«0Me» 


*iir* pW— I 

rl« »— * rl 


The doubled solutions appearing everywhere except at 
are necessary because from the factor ewwi or nnmi always combined 
with R^, through the factor or «'ni occurring together with 
co$ecH or dgu , simultaneous functions of iw+i)! and of (m— 1)2 
always appear, which necessarily require spherical functions with 
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the corresponding lower indices (m+1) and (tn-1) as factors. 

It has already been pointed out that the coefficients of the 
series pertaining to xmi« and are subject to certain con- 

ditions of a purely analytical nature, which is not the case for 
those of the series for g,H,z. Now there is another condition 
due to physical considerations, which naturally must be expressed 
in both representations; that is, that the integral taken over the 
entire earth’s surface of the force component perpendicular to 
this surface, must disappear. If we call the coefficients of 
SS^ieo$v)eoiml and JC («••)*»« mi in the spherical function series applicable 
for any function f of » and i, with 

Cl(/) and Slif)i 

then that condition is expressed in the simple equation: 

C*(Z) mm 0 

If we introduce the components r, jy, z , then by substitution of 
the expression specified for Z, we immediately obtain the slightly 
more complic<awed equation 

C*(5) + ^,(fl)+C;(^ -0. 

All these comments apply essentially when taking into afccount 
the flattening of the earth. The formulas to be applied in this 
case (which can be easily derived from the foregoing by introduc- 
tion of • ) , thus are: 

ayS "m -‘aX$inv.dyvcogi—07iinv.f^ee$tevmHl—‘TZ.tmvml 
ttfU -uX$inv.<tgvtinl-¥fi7nnv.f*&mcvcotl—fZ,9mvmi^l 

afZ mm 0-^.uXtmv — fZ.0mv 

aXrinv mm —^0eMv(S.tinveotl+H.tm9rinl)+Z.(l—mt^ 

0T$iHv ^ 0(’~S.$mviinl+H.im9eotX) 

fZ « —0r*{g.»invml’¥H.$mvtiitXi—Z.cetv 

The conditional equation to be fulfilled is thus: 

C’(^+^!W+ACJ(2)-a 

Since 4(=ryT+7*) is a constant, then the derivative of uXrinv^ 07tmv^ Z 
from S, H, z i.5 no different and no more complicated than that of a 
sphere. The converse, but practically unimportant problem, under- 
goes an important modification inasmuch as a closed expansion does 
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not result for s, H, Z ^ but for the products of these quantities 
with that is, The elimination of this factor is 

of course possible, but generally leads to infinite series. 

However, the factor occurring in connection with $Ttm» ^ which 
can be used in the form (l— causes only a small expansion 
of the calculation, without changing anything about its nature. 

Set-up and General Solution to the Normal Equations 

Everything is now ready to derive the normal equations. If 
fjjj . is one of the quantities belonging to the value v^ : 

mt1cm,innvi\ otEw^itinVi \ \ • Ti^i 

and if PSi (mit n ~ m, M+1, III+2 — )denotes the corresponding coefficients: 

A; 01; Dl: £l;il; &1, 

then the system of error equations runs: 

A* — JE^-Pl.J2l(eaM0 t — 1, 2, 8....2S 

For reasons presented in detail in B, p. 22/24, I have given 
equal weight to all these equations. In the case (not initially 
treated here) that no secondary conditions are to be met, we then 
have the following normal equations: 

(«»•«<) - ZK. J» — ■». »+l «+r 

or, in standard, abbreviated notation: 

[/-fia - 

nmm 

Here,v depends on the expansion of the series of spherical func- 
tions. For reasons given below, in the series for oJTnnv, the 
expansion shall be carried one step farther, i.e. »' is to be set 
greater by 1 than in the series for fiTsinv and rB. 

As mentioned earlier (p. 6), the system of normal equations now 
breaks down into two completely separate systems, one of which 

contains only F2, FZ** > and the other contains only FZ** 

as unknowns, because in general: 

SZ(tMVn-i) — (— 

and since consequently the sum [iCJBy, i.e. disappears for uneven 
values of (n-p) through an easily understoo<J abbreviation: 
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+ JC.U 

Let us write the two groups of equations in the form: 

Zrt*’'t«r’'J5r”] -t/«sr*i -»r*‘ 

z n**^*^. [«:♦*'*♦' s2**’*n - [/- k**’*'] - iT**** * - o, i. 2 . . . . 

or in shortened form, omitting the lower index m: 

(4) “ ?»■; “ f»«« *"0,1,2.... 

0 

For the numerical calculation on which tables I and V are 
based, we naturally set: 

?»■ “ Z (/■•,<+/»»,*-<) f»»+i ■“ Z 

Ot <•! 

The coefficients 'a' do not depend on the observation data, 
but only on the selection of the parallel circle to which these 
data relate. Thus, they can be used for every other calculation 
which is based on the same parallel circle. Naturally, its appli- 
cation is by no means limited to geomagnetic problems. Primarily 
for this reason, I have computed the solution of the individual 
equation systems, even though a frequent application of the obtained 
formulas (presented in the following pages) seems unlikely. In 
general, for similar problems, the Neumann method or a graphic 
derivation will be preferred (see A, p. 25, 26, and due to the 
reasons which Induced me to select this method, B, p. 21, 22). 

Regarding the coefficients ’a' presented in the following table, 
it should be noted that they have been derived from the values of 

computed initially by me and rounded to 8 decimal places, 
through multiplication with The reason for this is that 

after some delay and after a considerable part of the numerical 
calculations had been completed, I decided on the deviation from 
the usual method in the introduction of the functions R. I mention 
this because in a direct calculation of 'a* from R, the last 
decimal places were found not always to agree with those given here. 
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The differences are practically meaningless, which is why I omitted 
the time-consuming re-calculation. 


Coefficients of the Normal Equations (A) 


25.00000 


an — 23J2561 


Am • 51.67489 


m a 0 



- 1J7208 

an ■■ 

- 7.14895 

Cm * 

- 4.52886 

«M ■■ 

17.66068 

a«t » 

- 9.58886 

a« " 

- 8.44497 



a*t " 

18J8245 

«44 — 

— 4,WMff8 



• 

• 

a** • 

22.755M 

«ll — 

- 7.87997 


-10.15889 


- 8.75241 


16.18116 

«ii ““ 

- 7.51818 

«n ■■ 

- 5.58846 



au — 

21.88857 

an — 

- 8.12799 





•n ■■ 

21.95444 


M 

« 1 




Oh — 

4J5221 

Cm — 

- 4.55287 

an — 

- 6.18189 

an — 

52.78554 

a*t * 

- 7.64619 

a«i — 

-17.60242 



a« " 

88.15086 

a*4 • 

-15.60996 



♦ 


«M — 

89.05169 

a,t — 

56J6852 


0.89181 

at, — 

-11.89867 




44.68674 

«u “ 

-18.99684 





an — 

86.42590 


M 

— 2 




*•# “ 

50J5066 

a»» — 

5.59272 

a«t — 

- 0.47484 



an — 

58.58060 

««i " 

2.90698 





a« " 

49.65171 

an — 

56.80378 

att -■ 

6.22129 

a« — 

- 8 J5864 



«M — 

55.62099 

ati 

- 8.18960 





an « 

48.08946 


M 

i- 8 


• 


at* — 

49.27458 

an " 

4.69178 

an — 

0.59175 



an — 

58.52248 

a«« " 

7.04828 





an " 

57J5657 



an — 

54.98297 

at, — 

6.92004 





«n “ 

69.46148 


m 

» 4 






an — 

48.56728 

Cm ^ 

8.79217 





a« » 

56.99944 



an — 

68.27085 

«n “ 

6.16905 





an — 

59.45807 


m 

• 8 






a*t ■* 

48.09689 

an — 

8.11497 





a« “ 

55.48200 





•ii — 

52^)8705 

' 

m 

— 6 


On — 

47.76858 





•ii — 

51.18567 


m 

— 7 


On — 

47J2800 
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The expansions of the normal equations then run: 

■■ ^■Ia4>l,la4tfla4>l ^ 0| 1« 2 


The coefficients « appearing here are functions of the reported 
quantities 'a' formed in the knovm manner. Their logarithms 
rounded to 6 decimal places (which could be shortened even more 
for most applications) are found in the following table, as series shown 
for different limitations on the series expansion. 


Logarithms of the Coefficients of the Solutions of Normal Equations 

M « 0 


loffatt 

— 

9.007269 

logon 


9.088816 

lVo*% 

— 9.100947 

Ivon 

mm 

8.968866 




logon 

as 

9.466399 

IVOti 

— 9.882101 

Ivon 

mm 

9.261214 







Ivou 

« 9.483688 

logon 

•m 

9.227006 









Ivon 

rm 

9.209016 

logatt 

— 

8.704756 

logon 

— 

8.860648 

IVOn 

— 8.497846 







logon 

mm 

8.948616 

IVon 

— &741774 










logon 

— &980707 




legitH 


8.606616 

logon 

mm 

7.680816 









logon 

mm 

8.766960 






tcgUf 

— 

&602060 









Icgau 

•• 

9.642880 

logon 


9.700618 

hfgon 

— 9.627701 

Ivon 


9.420164 




logon 


9.819726 

Iv On 

— 9.716881 

IVon 

— 

9.616889 







Iv On 

9.678406 

Ivon 

— 

9.484078 









Ivon 

— 

9.827872 

logon 


9.065600 

logon 

— 

8.983011 

Ivon 

— 8.944128 







logon 

— 

9.191057 

IVou 

— 94)01887 










IVOu 

— 9.098688 




logon 

■B 

8.710228 

logon 

mm 

8.897748 









logon 

— 

8.869062 






l^on 


8.682167 



m w 

1 





log Ott 

— 

8.808738 

logon 

— 

6.766968 

IV Oit 

— 7.676061 

Ivon 

_ 

7.780408 




Ivon 

— 

8.419967 

IVott 

— a087296 

Ivon 

mt 

8.226187 







Ivou 

— 8.679862 

Ivou 

— 

8.831609 









Ivon 


8.629620 

logott 


8.294080 

logon 

— 

7.191096, 

Ugou 

— 7 J09080 






Ivon 

— 

8.293420 

Ivon 

— 7.674440 





logau 8.4SM88 


Lines 9-12, 19-22: The values of and %«n belonging to 

m=0 (not used in further calculations and presented here only for 
the sake of completeness) are somewhat inaccurate because they 
were computed from provisional normal equations whose coefficients 
«t* and *11 ■ deviate sometimes by several units in the last 
place, from the final values presented on p. 16. 


(A)* 
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log mm 8.290661 

lcy«M — 8J867tt 

kt9«H ■" 74168846b 
logan — a281828 


2cy«it — 8.278688 
kf «ii — 8.244891 

logatt — 7J06811 
l^mn<m a409722 

Iggan — 7.848844 
l^an 8J082482 
* aB2B641 

k|r«»i — ai87280b 

kf «ii « 8A44865 

Itgmn — a880884 



m -i 8 


log an aS02727 

logan — 7.286896, 
l^an 8.288679 

lagan — 6.486818 
kyaai — 7.014794. 

kyaaa -• a806808 

Ugan « a802627 

logan — 7J82867b 
logan ■“ at87ttl 


Ugan — a297998 



lagan 8.258166 

loga$t « 74188288. 

kyaa, -« 7.164920 

lagan « 8.260987 

log an — a261629 

l^an» 7.071821 
kyaai — aS69694 

lega,i — 7.299679. 
logan — 8.860114 


lagan — a248585 

■ 

« — 8 


lagan 8.810706 

logan 7.214404. 

kyaaa - a9916. 


ky«n — 8.242418 

ky «4i — 7J88176. 
kya«« — a847888 

log an ^ a810706 

k9>n — 7.814716. 
kyats 8.886006 


logan — a807878 



lagan — a266581 

logan ■> 7.882461. 
kyau — 8.282179 


logan — 8.260167 

■1 ■■ 4 


logan — 8.816919 

logan ^ 7.188986. 
7iyast — a846891 


logan"* a8l8667 



log an — a278761 

logan — 7.294812. 

kyaas — 8.2810tt 


logatx » a278610 



In several cases, conditional equations have to be taken into 
account. The resultant changes in the solutions will be given now. 

In the Z-series (as I will call the series used for the 
expansion of ) , only the condition — o has to be met. A modifi- 
cation of the computation only occurs here for the first of the two 
equation systems characterized by m = 0. The coefficients 
•Ml BMt . do not come into consideration or are set equal to zero; 
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the others (which I will call fi) obtain values whose logarithms 


are: 


l99$» — 9.188869 

Ugfn &897198 
h90n — 8.768808 


•I 0 

lot$n — aS49891 
U9$u « 9.060888 

— &816490 
IV/44 » 8.881470 


« 8.887790 
^ a747680 
Ug$m — &90668I 


The problem for the two equation systems belonging to m = 0 
Is only a little more complicated for the expansion of uX$inv and 
tTtinv . In this case, If F Is again to be set In series for 
B, C, D, E, In all 4 cases the conditional equations (B, p. 11) 


apply: 

where 


-0 

-0 




K2«+l 


denotes the value of the function r” at the North pole (l.e. for 
«> -" 0 ) . The coefficients of the sought solutions computed with 
respect to these conditional equations will be called r. Their 


logarithms are: 

M 

« 0 


logrn — 8.665270 

lopru - 7.829511, 

lop rot M 6.878826, 

— 7.711288, 


loprn &493678 

l^ru — 7.928064, 

» a017281. 


lopru ^ 7.964750, 

l^r** — a423535 

— a220089, 
lopru a886278 

a640271 

lopru^^ 7.671407, 


• 

loprn * &417877 

hpfu — 8.216791, 
lopru ~ ai4094i 


topr- * a619842 

lopru « 8a70867, 
l^rn — 7.820878 

4 


lofru"‘ 8.514684 

Uprix — 7.516348, 

lopr%\ ■“ 7X197879, 

lopru » a058681. 


lopru “ 8.519360 

lopr%% ■“ 7.965791, 

it^ ru ai20704. 



lopr%% — a416568 

4^ ai89689, 

lopru 8486642 

lop fit — 8.447558 

lopni - 7.941971, 

{<yri, - 7.88408(IW 



l^rtt &426269 

lopru -> a228196U 
lopru — 8.2889M 

• 

lopru — aS91718 

lopru - 8.207724, 
l^r** — 8X128786 


• 

lopru “ “O® 
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The solution for the X and Y-serles is much more complex, 
provided «>0 , because in the conditional equations appearing 
here, the coefficients of these two series are not separate. 

If we state in general: 

- 

with .• -.-«*!(•+•«)* j 

then the conditional equations are (see B. p. 11): 

‘A,-*!;,-©, •q.+'D.-o. *c;,+«D«-a 


The importance of al is that ai^tinv* represents the value of 
Si, for infinitely small values of *. The logarithms of the ^ 


coming into consideration here are: 


h0*\ — 0.288561 

loycjf ■■ 0J116SS 


— 0.287016 

i- 1JW2696 


* 0J20488 

— 1.142215 


— 0.8M065 

Icy*; — 1.252522 


* 0.588046 

lcy«« • 0.977121 


logm\ » 0.7095M 

— 1.2M^1 


» a797610 

l«f •• IB 1417M1 


\ogm\ . 0.866762 

mm 1.665782 



U0674S 

L417061 

L6470K 

L2180ei 

13767et 


— 1J1108 


The simplest method would be to determine the coefficients of 
both expansions linked by a conditional equation, by means of a 
common balancing. Accordingly, I have done this by suggesting a 
detour which permits a certain judgement about the reliability of 
the final results (see B, p. 36, 37). The cited conditional 
equations are of a purely analytical nature; they tell us that the 
horizontal force at both poles is unequivocally determined in its 
magnitude and direction. They would have to be inherently ful- 
filled if the coefficients of the series for X and Y were computed 
on the basis of our knowledge of the force distribution over the 
entire earth's surface and if this distribution were expressed 
without remainder. Now this knowledge is missing for the calcula- 
tion to be performed here, for the two polar spherical indentations 
(the other side of 60° N. and S. latitude) and related to this. 
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the computed, first coefficients of the series expansion are not 
dependent on the others, which are neglected. Consequently, 
the numbers obtained In the Independent calculation of the 
coefficients of both series need not necessarily satisfy the 
conditional equations, and they will also In general not satisfy 
them. The amount of the remaining error apparently permits a 
view of the level of reliability of the results. 


Consequently, I first computed the values B,C,D,E Jndependently 
and used the general solutions with the found coefficients «. 

To the found values (printed in B, p. 54, 55 under II), which I 
will call B’ , C, D’ , E', I have added those corrections (after 
calculation of the remaining error of the conditional equations) 
which make these errors disappear and the sum of the error-squares 
of the original equations are brought to a minimvMn. (The final 
results thus agree completely with those which would result from 
a direct, joint balancing using the least squares method — this 
must be stated in order to prevent any misunderstanding). 


For the sake of brevity, I shall be satisfied with a statement 
of the numerical results without their somewhat cumbersome, but 
easy derivation. 


By insertion of the corresponding values of B' and E', let: 

- A. , • Ji 

Then we have, depending on the expansion of the series, the follow- 
ing equations: 


. ^/-[7.5448S]A, 
S\ ^'^[8.05062] A| 
Si — ^'-[8.22918] A, 
SI * ^'-[8J25i5lAt 

S\ m. Si'^[7.780a7] A| 
S; — ^'-[8.81780} A| 
St - J*'-{8J6646]Ai 

Si mm. £]'~[8J0608} Ai 

Si - 


- £*'+[7.8407:341 
£» .. r*'+[7J814QAi 
£* — £*‘+[8.088211^41 


£J - £f+p.9027434, 
£•- £*'+[8.4068734, 


£• - £^+[8.787183 4, 


j* .• 

^ i. 

J* - 

Bi -> ^*-[8.8M12j^, 
B! - ^-[8J8«S;|jrt 

B! - 

J; * J^-t7.W7M)4, 
JJ - 5S'-[7J4«»]A, 
J* - 

Bl - ^'-t 7 .M 70 SlA| 
JS - ^-[ 8 JMf 734 f 

Bt - ^-[7J2l4«]£b 
J* aj'-[7.7»1473li 
Bl — £*'-[ 8 J 1269 ]J& 

Bl - ^'-[7.«068«3ft- 

- ^'-[ 1 . 12296 ] A 

SI — £^-[ 2 J 26 »l 3 J& 

. Bl'-[MmO]A» 
Bl - [7.48269341 

Bl - .B2'-[81099634t 

» j;*- [7.696793 A, 
£* . J**-[&767103it 

B* - [7.600883 £s 

£• - [8.460873 £• 

Bl • .B;*- [9.181673 JT, 

- SJ'-[ 7 J 05263 a« 
SI - J;'-[ 8 . 671803 a, 

j; - JB;*-[7J70863£« 
Bl — J*'-[8J29683£« 

flj - JJ*-[9.09196J£< 


i 8 J - 

El - JJ'+pJ 89110 )^ 

z; - lj'+CM*aw 3 ii 

El . JE{'4-pJ818q« 

i; - jr+[UMtQCi 

El « J!^-[94n8»3« 

j; . £^4>p^Mie23At 
^ . £»'+[ 8 J 8 & 4434 , 

El » Jer 4 -[M 108 ^ 6 , 


El - ir+[«^«3*i 
£J - £** 4 P^ 069 }« 
- jy+pJiOOl] Jb 

JE; • £rH-p48i883« 
J!« - £S’+[M7778] Jb 

El - £T+[M>M03« 

£» « £r+pJOW 83 A, 
El - £S'+pJ 618 q 4 , 


r, - JEr+[MB 29 q 4 . 


£! - £;'<f pJ 28123 « 
El - 

js; - £r‘(*[i>^ 763 A 
£• - jE:'+[ 8 J 178 i 3 A, 


£» «. £*'+[ 6.87114 £« 
£J - £r+[« 4 > 8774 £,. 

£> - £;'+[ 8 J 11413 £« 
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A glance at Che coiKilClonal equations teaches that the 
coefficients B and E can be replaced in the reported fornnilas by 
C and -D. This substitution naturally must al so be performed in 
Che quantities a and E, so Chat these can be defined by the 
equations: 
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Nun>erlcal Resolution of the Normal Equations to Derive the Series 
Expansion of mXtImvt /fF«iiieand r*- 

Now In order to obtain the coefficients of the series for the 
three force components, one only has to Insert the values In 
table V, which express the state of the geomagnetic force at the 
earth's surface, Into the formulas reported above. 


First determine the values of f. For any value of m, we have; 

where stands for the earlier cited 6 quantities 

etc., see p. 11). The ccxnputatlon provides the following numbers: 



ahtinv 

aKHnv 

0ltinp 

0L9mp 

rn 



S4S744 


> 1821 


6017 



-nsaeo 


6M 


84986 


If.m 

> 40641 


> 68 


— 89670 


wiia 

> 16164 


88 


- 1068 


i/.ipj 

10864 


1789 


6869a 


I/.ilS 

— 28077 


> 18a 


-611774 



14666 


> 607 


mynif 


If.Ki 

8664 







~ 66828 

19809 

>176276 

-49806 

16009 

-866690 

LT.SU 

86088 

-22669 

- 9967 

-84184 

116928 

- 1070 


9883 

186 

17928 

I88tt 

- 80688 

47667 


- 4118 

8416 






11829 

101860 

I87tt 

78567 

-217888 

67884 


- 44449 

10278 

> 8276 

8967 

— 400 

- 16878 


20798 

>89197 

> 8668 

-17889 

68694 

- 18069 


- 87484 

- 8444 

-62197 

26888 

- 68887 

- 96111 

i/.iia 

18907 

2614 

- 1611 

84107 

- 81802 

960 

C/.i^ 

18467 

- 888 

- 1996 

- 886 

6479 

- 1801 

Lr,aa 

- 7892 

88661 

1266 

618H 

-187010 

- 8466 


80810 

- 2684 

28 

19640 

- 86416 

4917 

CA.i^ 

44 

694 





i/.«a 

8268 

7948 

- 86870 

26864 

- 28841 

- 640M 


- 8748 

- 8841 

- 206 

16a 

- 4464 

- 9907 

l/,S}} 

- 188 

> 866 





ir,sQ 

7478 

18899 

114M 

-188M 

81899 

14917 


8796. 

1898 

6871 

- 6686 

28741 

18806 

C/.JQ 

8660 

8819 

12^ 

6071 

- 6669 

6801 

C4J?a 

- 8808 

116 

- 8077 

4478 

7881 

11900 

ir.sn 

6460 

> 961 

6900 

-4816 

- 4766 

6171 


- 49 

- 176 
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By substitution of these values into the formulas of the 
preceeding section^ we obtain the desired coefficients of the series 
of spherical functions needed for representation of Utinr^fE, 

Their values are found in table VI for three different limits of 
these series. A comparison of the values of the first, comiron 
coefficients in these three cases permits a certain view of the 
attained approximation of the true values of these coefficients. 

The occurring differences are a result of the fact that the 
valuable observation material does not include the entire earth's 
surface. If this were the case, then the individual normal equa- 
tions would contain only one of the coefficients and these would 
then not be interdependent; the computed ones could thus not be 
affected by the ones neglected in the expansion, and change with al- 
ternating extension of the expansion. It would thus also be 
formally possible to reduce the mentioned differences by taking 
into account the unused observations in the polar regions. As 
long as our knowledge of the geomagnetic force distribution in 
these regions remains so deficient, particularly in the Southern 
hemisphere as is presently the case, only little more than an 
apparent increase in accuracy would be attained by this. 

In order to have a completely well-defined analytical expression 
for further calculations, I rounded off the computed coefficients 
where possible, to whole units of the introduced unit </» . This 
is the case without exception for those of yZ. For those belong- 
ing to X and Y which are linked by a number of conditional equa- 
tions with irrational factors, one of those appearing in the same 
equation as a fwiction of the others (all of which were rounded off) 
would always have to be ccmiputed aikt thus would not permit any 
arbitrary modification. These coefficients are given in the table 
to two decimal places. Originally I had selected those of the 
lowest order (see B, p. 23, 26). The values computed under this 
assumption form table VIII of my previous report (B, p. 57). 

Later it seemed better to choose one of the latter coefficients 
of each conditional equation belonging to the highest order, as 
a function of the others. On the one hand, the corrections occurring 



in Che rounding are smaller in this case Chan for Che flrsC meChod 
(as a glance aC Che condicional equacions will show); on Che 
oCher hand, Che coefflclencs of Che flrsC orders are chose which 
analyClcally define Che represenCed sCace, whereas Che very flrsC 
one cannoc he given accuraCely as a funcClon of Che following one 
and ic Is of no Imporcance for Che deflnlclon of Che sCaCus. 

Thus, Che coefflclenCs of Che X and Y-serles compiled here 

In cable VI^ do devlaCe In pare from chose of Che earlier reporc 

(B), and Che same applies for Che compuCed coef flclettCs of Che 

series for U, W, V and 1. (However, due Co an unClmely discovered 

failure In B, Cable XIa, b, p. 60, 61, chose values of Che 

c<^puCed coefflclencs k^...K^, 1^...L^ were prlnCed which resulC 

from Che quanclcles B, C, D, E rounded off by Che second method 

and reported here. Only for k and 1 are the coefficients B and 

o o 

D based on the first type. Let us take this opportunity to 
correct another minor error which showed up in a repeated, 
thorough examination of all figures in Cables XIa, b, c. The 
values of belonging Co r m su* and r -= lou* should be changed Co 
-11000.7 and -12784.3). 

In Che following Cables Vila, b, c we now find Che values of 
the coefficients of the trigonometric series compuCed from the 
numbers of VIj which represent X, Y, Z, for u * 0°, 5°, 10°, 15°.. . 
180°. The k^, 1^, m^ independent of the geographic location, are 
given directly; all others whose numerical values are round In B, 
are given by their logarithms. They were computed to hundredths 
of r and Chen rounded off Co tenths of this unit. 

The Cables Villa, b, c based on this, finally give Che values 
of Che force components X, Y, Z In whole units r for all points 
separated by 5° longitude (in Che calculation, tenths of a unit 
were carried over). 

The three Cables VIj^, Vila, b, c. Villa, b, c contain the main 
resulc of Che present Investigation. They present Che same dls- 
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tribution of geomagnetic forces in three different forms. How- 
ever, only VI^ (which contains the coefficients of a series of 
spherical functions) defines the distribution without additional 
information, because accordingly the calculation can be carried 
out for each point of the earth's surface. The numbers of VII, 
the coefficients of the trigonometric series belonging to 35 
parallel circles (except the poles), define the state of 
the force initially only for all points of these parallel circles 
(and for the poles). In VIII finally, this state is shown for 
2522 separate, regularly distributed points. All three represen- 
tations are theoretically entirely equivalent if we add the con- 
dition to VII and VIII that the force distribution be expressed 
by a 6th order (for X, 7th order) series terminated with spherical 
functions, and for VIII the added condition that the trigonometric 
development along the geographic longitude be terminated with 
the functions of the 4-fold angle. 

The state defined in three ways is also characterized by the 
fact that of all possible states, it most closely approximates 
the state observed by the numbtis illustrated in tables III 
(or IV) — a statement that has a certain, though practically in- 
significant incongruity that it depends on a specific setting 
of the weightings. 

The main value of the figures reported here lies in the fact 
that they form a convenient starting point for any future compu- 
tation of potential based on new material. In order to perform 
such a computation, one will calculate that value of the measured 
element for each location where a valuable observation is made; 
said element results from the analytical representation provided 
here, and then the difference between observation and calculation — 
which naturally also contains the secu' ar revision — will be 
selected as the basis for a refined calculation. (See the dis- 
cussion by E. Sobering in Geogr. Jahrbuch, XV, 1891, p. 143/146, 
which is still accurate if we proceed in a non-Gaussian manner by 
using graphic mtchods or Neumann's formulas). 
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In order to determine the computed values of X, Y, Z (and 
thus also those of S,i,i) for the individual, randomly distributed 
observation points, it would be best to proceed so that first the 
tables VIII are expanded for each full degree of longi’-jde and 
latitude by interpolation, and then to pass to the observation 
point through a second interpolation. If no greater accuracy is 
needed for the computed value than about 5 to 10 r — which should 
be sufficient with regard to the accuracy of most observations 
in general — then in the second operation one can set by everywhere 
with linear interpolations; if more accurate values are desired 
for individual locations where particularly accurate observations 
were made, then under all circurstances it is sufficient to take 
the second differences into account. The first operation, the 
expansion of table VIII to degree-intervals, does require an 
interpolation with fourth, and sometimes even with fifth differences 
if the values are to be obtained to approximately 1 / accuracy. 

This rather complicated method for a table with double entries 
will be difficult to apply to the interpolation of a single value 
and is thus subiect to a significant simplification such that the 
computation is always to be performed for the same, very con- 
venient interval of + 1/5, + 2/5. The pertinent formula will 
be derived below and a practical example of its application will 
be given. 


Let A, be the value of a component belonging to the tabular value 
of u, let At and a, be the 2nd and 4th difference of the values 
standing in the same line , a',, a; and AJ, aJ are those of the 
1st, 3rd and 5th differences in the preceeding and following 
spaces. Then we obtain the values belonging to (m— 2 “), («— i*), (m+2*) 

of the quantity to be interpolated by substituting in the expression 

A J-JL I , I ^ a;+a; 

1 ■ 2 1.2.8 ’ 2 1.2.3.4 1.2.3. 4. 6 ‘ 2 


for n, the series using the values -2/5, -1/5, +1/5, +2/5. We 
thus obtain the four functional values: 

.4— (J'+jB*), c— (iX+i)'), e+ijy+D-), 
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if we use as an abbreviation: 

^ " T^*“” 250 ^*^1556^ 

^ “ ^■**50^“ " To^‘“i26^‘*’imo^ 

(The factors of 5« and a, can be rounded off without causing any 
notable error if one does not wish to prepare a few small 
secondary tables. Equivalences arer^gjgt jjig'gisw i56’ fio* 82o‘) 


As an arbitrary example, let us use the calculation of Z 
for a few points of the 15° East meridian. The difference out- 
line and the compilation of the auxiliary quantities obtained 
from the differences then look as follows: 

■ ^ A. ^ ^ ^ Ttr-t-r) i>_ 

20* 50A06 

-1771 


25 * 48686 — 12 

—1788 —184 


80* 

46852 

—196 

-34 





—1979 

—218 


50 

—889.4 

—1944 

85* 

44873 

—414 

16 


448394 ± 862.1 

44864.7 ±4804 


-2393 

—202 


43 

-472.7 

-236.0 

40* 

42480 

—616 

59 


424304 ±10704 

42467.8 ±534.1 


-3009 

-143 


41 

—5974 

— 2M.8 

45* 

39471 

-759 

100 





—3768 

— 48 





60* 

857(» 

—802 






-4570 

56* 81138 

From this, we have for u = 33°, 34°...42°, the following values of Z: 



S3* 44839.8+862.1 — 45702 38“ 42430.4+1070J — 48601 

34* 44864.7+430.8 — 45295 89* 42467.8+ 53U — 43002 

35* 44873 40* 42480 

36“ 44864.7—430.1 44435 41* 42467.8 — 5344 ■■ 41983 

37“ 44839.8—862.1 ■* 43978 42* 42430.4—10704 — 41860 

To check the calculation it is simplest to use the difference 
series of the found number series. In the present case it turns 
out that the second differences — which alone have any significant 
value — are sufficiently regular. The small anomalies in their 
profile can be attributed to rounding errors. (It may be important 
to note that the irregularities caused in some cases in the fifth 
differences of the original series, are intensified somewhat by 
a special circumstance. Strictly speaking, the numbers cited in 
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the tables are not exactly equidistant because they (see p. 4) 
are derived from the values of v rounded to whole arc-seconds; 
their attendant u differs from the round numbers used as arguments 
for the tables. This difference Is a few tenths of an arc-second. 
Practically speaking, this inaccuracy Is meaningless; for none of 
the computed values of X, Y or Z does It cause an error of 0.5r). 


Perhaps even more convenient, though It contains three 
successive interpolations, is another method which Is somewhat 
better than the above method, at least for regions with numerous 
and accurate observations. It consists in first finding the 
functional values for the middle of the 5°-intervals by interpo- 
lation, and then interpolating for every The first operation 

where the differences of uneven ordering numbers drop out, need 
only be carried out to the fourth difference, the second operation 
where the above-developed formulas are used, is carried out only 
to the third difference. The last interpolation within the 
interval which leads to the values for the individual observation 
points, is performed linearly only. Through the repeated inter- 
polations (longitude and latitude) a pile-up of rounding errors 
occurs, in addition to the errors due to neglect of the higher 
differences; the total uncertainty of the final values will then 
not generally exceed 2 to 3 r • 


If A, and a; are two sequential functional values (belonging to w, A 
and n+6*. A or to it. A and m, A+S*) and a,, a;, a,, a; are the second and 
fourth differences standing in the same line with them, then 


the interpolation for the middle of the interval gives the value: 

a,^.a; 1 a,±a; A4-FA1 

S 8 2 ^128 2 


or when 





r 


is used , ^(P + P' ) . 


The next interpolation occurs, as mentioned, by using the 
earlier formulas which are simplified by elimination of a,. We 
thus have: 
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A ■■ I S 




A,* 


J> 



and the following values must be formed: 

C^(iy+D'), C+(2)*+i)'), A+(B^+B*) 


An example may serve to illustrate the method. The values 
of Z reported above give the following values of P for u = 30°, 
35°, 40° and 45°: 

46876.7 44926.1 42668.4 89668.4. 


The average values of the sequential ^(P+P') are: 

45900 43742 41063 

and denote the values of Z for 32.5°, 37.5°, 42.5°. 

The following outline contains an additional interpolation 
for all intermediate points from to ^°. 


• 







80% 

46862 








— 952 





82% 

46900 

— 76 







-1027 -29 


-204.4 


—102.4 

35% 

44873 

-104 

44864.7 

±480.0 

448704 

±2164 



—1181 -27 


-226.4 


—112.7 

37% 

48742 

-181 

43731.t 

±477.1 

487394 

±2384 



—1262 —24 


-261.7 


—126.4 

40% 

42480 

—166 

42467.4 

±684.4 

424764 

±2874 



-1417 —20 


-2824 


— 14U 

42% 

41068 

-17fr 







-1692 




. 

46% 

89471 







The following table results from this; it fully agrees with 
the results of the earlier calculation: 


u Z 

84% 45296 
84% 46086 
88% 44878 
86% 44666 
86% 44486 


t$ Z 

86% 44209 
87% 48978 
87%. 48742 
88% 48601 
88% 48254 


% Z 

m 48002 
89% 42744 
40% 42480 
40% 42210 
41% 41988 


The differences formed to check the calculation exhibit a 
quite regular profile. 
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Calculation of the Coefficients of the Potential 

From the series computed above for and /Ftmv , the 

potential of the horizontal force, provided one exists, can be 
determined and illustrated likewise in a closed form. For this 
purpose we have to calculate the functions: 

* I 

£7 — J aXdv E^+/(r,i), IT — //Fwiirtfl ■■ 

t « 

which generally each contain a part (/(r,i) and which cannot 

be represented by a finite series of spherical functions. 
denotes the i -independent part of which can be expressed by 
spherical functions (see B, p. 9). 

If it turns out that U = W, in which case / and * always 
disappear, then the entire magnetic horizontal force at the earth's 
surface can be defined by a potential which is determined by 

F- bU^bW 

g 

with b = 6,365 * 10 cm as the polar radius of the earth. In fact, 
we then have: 

y _ 1 8W ^ 1 iV 

« “ *5 »r ’ “ /ttinv 9JL “ /iFlmp IT 

as it must be according to the definition of the potential (see A,p. 7). 

If U and W are not equal, then that part of the force to which 

a potential is ascribed, remains undefined to a certain extent (see 
A, p. 17). If we take it to be as large as possible — which is 
evidently useful, even though not entirely sufficient for an 
unequivocal definition — then in the simplest case we set: 

F - |(J7,+ irj 

To characterize that part f the horizontal force to which no 
potential will correspond, we use the statement of the difference 
(W-U), through which ^Wq~Uq) is given with consideration to the 
'a priori' specified form of / and z. 

The determination of U now takes place by means of the follow- 
ing formulas (see A, p. 20; B, p. 12). 
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As abbreviation, we write: 

aZiinV mm Z-liiSi* i.e.^i 


and let: 


Jffiteetml+CttHnml 


so that: 


nt 

n. 


xL - 1, 


A ‘9 ^ fin+J>) (2||H>^1) (^1) 

■ *" <•»+!>- 1M2im+2j>— 8) (2m+2j»—l) 





— V 

1 ^ 1 • • 

■■ •g'*’— -j- «m8ff 

" + g «n 4 p 


n% ■« 1— pMp 

8 8 1 
^4 ■■ y — -j* «*p+g «*8 p 

» 8 8 8 1 
“ 15 — 1* «»w+^ eo*8p— gg «Bf 8v 


Then: - Z»«i7«+ZPiai 

if we compute the quantities: 

rtm » faCMIR^ + CafinMa, OM «2 + -HI «» mA 

from the following equations: 

;r« - #iI-42+#.Lir*+#.i.Jr*+ .... 

4/s^' - .... 

pIPT*- 


pLp 2 -A-c^'+/4ir*+A-<r‘+.... 

kLps^»- #4i.4:^*+#.i.-iz^‘+ 

pLp:^- ....' 


For m = 0 it is easy to see that n% equals the corresponding 

constant «? introduced above (p. 14) which appears in the conditional 

equations for the coefficients and c”. By virtue of these 

o o 

conditional equations, we evidently obtain: 

IT, * 0 , P{ « 0. 

The other equations run as follows after introduction of the 
numerical values of » and p. 
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F\ — [0.2870156] J2 + [0<4146S18] iLS + [0.4M602S] its 
Fi — [9.8204880] -f [9.9008896] + 

PS— [9.6006468] IS 4- 

JFS — [9.9988827] A\ + [0.0920806] ^S + [0.1698794] Al + 
FI — [9.6967881] + [9.7641869] A\ -f 

[9.6217680] .iS 4 


« [0.2385607] A\ 4 [0.8856660] 4 [0.4786588] A\ 4 [a6448627].4S 4 • • • • 

F{ » [0.0194590] A\ 4 [0.1124568] A\ 4 [0.1781667] 4 • • • • 

Ft — [9.7056520] A* 4 [9.7718614] ij 4 . • • • 

— [9.5282454] 4 •••• 

FI ^ [0.8494850] Af + [0.4655598] A* 4 [0.5426626] .ij 4 • • • • 

Ft « [9.8845088] A* 4 [9.9116066] .4* 4 • • • • 

— [9.6066640] ^• 4 ... 


», — [0.2870166] A* 4 [0.8494850] .4S + [0.4170618] .4* 4 • • • 
Ft — [9.8829588] 4* + [9.9505246] 4 • • • 

PS— [9,6261229]4S4-. 

P| — [0.1216190]4S4[0.1950578] 4S 4 [0.2558220] 4S4--- 
PS — [9.7846480] 4S 4 [9.7949127] 4] 4 • • • 

if— [9.6402698] 4S 4 .- 

nr, -t [0.3204890] 4* 4 [0.8860854] 4S 4 [04849888] 4] 4 • • • 
PS — [9.7986976] 4S 4 [9.8426465] 4S 4 ■ • 

Pi— [9.5658860]4|4-- 

PS — [0.0000000] 4S 4 [0.0454098] 4S 4 .... 

PS — [9.6680161] 4* 4 


— [0.8460658] 4* 4 [04817082] 4S 4 ...^ 
PS— [9.7281829] 4S 4.* • 

PS — [9.9186864] A* 4 [9.9441904] 4 ] 4 ■ - • 
PS — [9.6076091] 4J 4 . • . . 

[ See , note below>*f ] 


*The numerical formulas compiled here are not given in the form 
used in the general presentation in order to simplify the usual 
logarithmic computation. But they can naturally be converted 
into this form immediately and this would be an advantage for a 
numerical calculation when using addition logarithms. 


33 


ORIGINAL PAGE « 
OF POOR OUALITY 


These equations whose coefficients were given as more highly 
accurate, permanently valid values than would be necessary for 
present purposes, show that the series for aX$in0 has to be 
developed down to 7th order terms if those for U are to be 
obtained to the 6th order (equal to those for Y, Z, W). 


Now if we set the numerical values reported in table VI into 

the above, general formulas for the coefficients a” (i.e. 

m 

Clirinmi) > then we obtain the coefficients of U specified in IX. 

For the case of the farthest expansion of the series, we also find 

these values printed in B, p. 58, table IX. The small differences 

existing in the last decimal place between the two statements are 

due to the different rounding of b’^ and (see p. 18). 

m m 


The calculation of the coefficients of Tf— from those of 
07tinv is very simple — in contrast to the derivation of U above. 
Evidently, since: 


and 


then: 


IT— ““ — 


■ «l «■! 


The function ^(e) which appears as an integration constant 
with respect to i, is as already remarked, the part in free 
of A: 

^(f) — ziF* 

••I 

By means of these expressions, W breaks down, like U, into a 
finite series of spherical functions and into a part which 
cannot be represented in this form: 

*(t).A--Ar'c:j: 

These latter, as well as the coefficients of W^, are also 
found in table IX, which further contains the coefficients of 
^.(Uq + Wq), i.e. of V:b. 


ORlGiNAt PAQEJS 


The two functions U and W or the equivalent expressions: 

F-ir- (F,-rj + x{r).a-/(,,i) 

together determine the state of the magnetic field in the earth’s 
surface uniquely and completely, as is done by aXftnv and fiY$inv 
together. 


V, the potential of the horizontal forces at the earth's 
surface, can be broken down into two parts, V. and V , by means 
of the expression found for ^Z\ these parts represent potentials 
originating from agents in the interior of the earth and those 
in outer space. This is done by means of the formulas (see A, 
p. 23; B, p. 13): 


7 — + 

■■ ft (rl cot Ml mail} XI 

with ^ A " 

tm * A “ 


The constants #1 and ^ appearing herein, depend on the 
flattening of the earth. For the case of a sphere, we have: 

" 2ii+i ^ iii+1 

The only slightly different values which result from using 
the flattening derived by Bessel (1:299.1528), are found in B, p. 
46, 47. Their logarithms are contained in the following compila- 
tion. 
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0 


0.0019S84 


9.5t4041l 

».5i5SOC4 


9JOS415C 

9.S0SC90S 

9.9035999 


9.1569999 

9.1564595 

9.1569999 

9.1574967 


9UM71997 

9.0479690 

9.0474949 

9.0479716 

9.0464010 


9.9600494 

6.9600967 

9.9609519 

9.9604954 

9.9609494 

9.9619914 


94975019 

94975979 

94976497 

94979199 

9.6990650 

9.9999999 

9.9997707 


0 

1 


•— OD 


94940411 

9.5999937 


9.6094799 

9.6099689 

9.6016499 


9.6999799 

9.6399169 

9.6390944 

9.6316994 


9.6490059 

9.6479790 

9.6476999 

9.6477591 

9.6475599 


9.6577969 

9.6577119 

9.6576707 

9.6575M6 

9.6574991 

9.6579589 


9.6649909 

9.6649999 

9.6649N1 

9.6649599 

9.6641909 

94641110 

9.6640199 


The resulting series derived for V. :b and V :b have the 
coefficients specified in table X. The same table also contains 
the presentation of a function afibi obtained from (W-U) and which 
conversely is sufficient for the determination of (W-U) so that the 
data contained in X provides a new, complete and unique represen- 
tation of the geomagnetic field. According to B, p. 13, we have: 

4/ifm V ip H ^ bnbtmv [ 44 


Jp 


The function i introduced here, means that an electrical 

current penetrating perpendicularly into the earth's surface 

having surface density i, would generate just that part of the 

magnetic, horizontal force (also expressed In (W-U)) which is not 

generated from the potential V. Since the unit of the numbers given 

4 -1 

in X is equal to o.i‘e»i-»^»#“« , i.e. 0.1 Ampere • cm and 

Q 

furthermore, b = 6.356* 10 cn, then we obtain from these figures 

the current strength (or rather the o/-fold multiple of it which 

2 

differs insignificantly) in the unit Ampereikm , if it is multiplied 
wlthaiMO**: (6.856. 10*), i.e. with 1:635.6 or 0.001573. Note also that 
positive values of i would indicate a downward-directed flow. 

To derive afibi from (W-U) we use the following formulas whose 
numerical coefficients which are again given to 7 places for the 
sake of overall accuracy, even though at the present at most 4 places 
are needed. 
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] dn* * 

— (OAMMona; 

1 dS* 

- -tO4ta»4i0]iJ-ia77*0B4(qJ! 

* dS^ 

klti if " JJ-CaWTilTtJJiJ 

1 dR* I 

1 ds* 

-{o.7»«»jS5-(at79MW]aS“P-«ni€iO)j; 

“ [(MNXKKNX^fjMlH 

1 dR* 

mni if "* 

J5T5 ^ " »•*?■•- lo.*MO*463 Ul -{atisesM] ^ 

1 <IJ?* 

■■ loJSiieoi) «#••!«»-* 

1 

•* [0JU«M7]«D»e#mir*-[0.7f761W]^ 

1 

“ I1.1W7420] emtmn r-» ~ [0.9978176] Rf - [14)066^] R* 

> 1 

;sr; IT “ U^I0«*47]-[0«96708]fi! 

^ - [1.5812910] - [1 J5M198] - [14)856780] .BJ 

4^^-[04890456]a». 

1 dR* 

H5-J? - |i M»:n»]»tr-tAii*imeJA5 

1 tfJ?* 

iS^lf " [1.7180618] (wr- [1.2815178] - [1.1082442] it* 

1 dn 

5^,-jJ - fAOOOOOOO]*. 

_ [i.jHai,jrt..-ii.o:»itii]ji! 

1 tfJ?* 

Pi?«»41725]mr- [1.6457700] .»;-[!. 140187l]i^ 

1 dS* 

tmvlf “ [0.7976108] «*iwr 
1 dS* 

•nii "aF “ «w * «*• r ~ [14»47276] 

*Change 0.2385607 to read 0.2385606 and change 0.8996702 to read 
0.8996703. 
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It is easy to see that when substituting these expressions 
into t the coefficients of cMpmnv^u ■ ■ mpmn^ 


become aggregates of J9S;, oi;, X)Sit -El ^^i^^ disappear due to the 

identical conditional equations valid for the expansion of 

mX$inv and 4Tmnv . For m0hi we thus have a finite series of spherical 

functions comprising the first five orders here. That the 

coefficient of is equal to zero therein (likewise as in the 

o 

series for V. :b and V :b) expresses the fact that the algebraic 
sum of all the currents penetrating the earth's surface, i.e. 
the integral of Hk taken over the entire surface, disappears. 

The surface element 4 w is equal to a 0 i* tin vdvdi , This explains 
that the expansion is obtained not for 1, but for «t, i.e. 

— except for constant factors. 


Due to the results reported above, the problem to be solved 
has been completed. However, we still have to investigate the 
validity of the obtained results. Since this has already been 
done in sufficient detail in my preliminary report (see B, p. 34-43), 
I do not believe it necessary to go into detail again here, but I 
will limit the discussion to a brief presentation of results. 

It turns out that the presentation of the force distribution 
defined by «x«tnv, /Ftmt, can be viewed as generally satisfactory, but 

that the derived functions V^:b, Vj^:b and adbi are affected by such 
considerable uncertainty that only the first, absolutely la-ge 
coefficient of V^:b can be viewed as determined with sufficient 
accuracy. and m$H however, which assume small values only, 

can hardly be viewed with enough certainty to believe in their 
existence at all. The reason for this considerable uncertainty lies 


mainly in the fact that the individual series coefficients cannot 
be calculated Independently of each other; this would require a 
knowledge of the force distribution over the entire earth's surface 
which we do not now have. Even though I still felt justified 
in giving a positive result in my preliminary report, this was 
done in the final analysis only through assuming very large and 
broad systematic errors in the observations of quantities to be 
explained and regular distribution of the value of 1^ (see B, p. 

43). Thus this decision is based primarily on the same principles 
used by L.A. Bauer ( Terr. Magn . Vol . II, p. 11) and v. Bezold 
. Berlin Sitz.-Ber . f. 1887, p. 414) in their discussions. The 
investigations since published by Schuster, Rucker et al . reinforce 
the weight of the counterarguments so much that the real exl stance 
of V and of 1 have to be designated as at least quite doubtful, 
and chat Che reliability of Che empirical fundamentals to be 
evaluated in any new calculation of potential appear in fact to 
be significantly affected by systematic errors. This applies in 
particular for i which can be determined Independently for each 
geomagnetical ly accurately investigated part of the earth's sur- 
face, and its calculation using the Rucker method now permits 
a check of the results obtained here (see A, p. 16). Just this 
possibility for a mutual verification — which would lose its 
significance only if we had a complete empirical knowledge* of 
Che distribution of the earth's magnetism — itidlcates the desirability 
of refining the expansions given here for the entire earth's surface; 

this would also be needed for a more accurate determination of V_. 

a 

The most important task of future research related to the 
determination of the spatial force distribution of the earth's 
magnetism, must deal with filling the large gaps, especially in 
the polar regions, and then on the oceans and continental interiors. 
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IVa. Coefficients of the Trigonometric Series for the North 
Component X 
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IVb. Coefficients of the Trigonometric Series for the East 
Component Y 
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IVc. Coefficients of the Trigonometric Series for the Vertical 
Component Z 
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Vb. Coefficients of the Trigonometric Series for 07mmv. 
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Vc. Coefficients of the Trigonometric Series for yZ 
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Vila. Numerical Values or Logarithms of the Computed Coefficients 
k and K in the Series Expansion of X 
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XT 

99 

—134 

38 

604 

—937 

8®S 

—479 

- 75 

786 

f — 43 * 

—••5 

*69 


XI ' 

S» 

7«7 

797 

1*3 

— 93 

-l8t 

— ja 

1 

—763 

1 — 4*5 

1 

— *74 

•331 

t 

1 

AJ 

— »35 

— 34 

S «4 

— 3*3 

—1*1 

35 

s6i 

-96 

-533 

1. 169 . 

34 

—404 

;<r ’ 

XT : 


«9 

-483 

7 ^ 

—539 

544 

- 544 ' 

*50 

•65 

; -597 

too 

— 161 

j 

XI 

- 79 * 

$®« 

-69* 


479 

^56 

537 

*34 ] 

1003 

1 7«7 

m 

Sio 

t 

XX ^ 

4*9 

>5 

— *80 

16 

-••7 

*76 

374 

—40* 1 

7 r 

i 4*3 

738 

-84 

9 T ! 

XT ! 

*07 

190 

-733 

709 

-433 

403 

-363 

141 1 

—318 

; - 6$9 

1*86 

-793 

1 

XI ■ 

-937 

1131 

-7*0 

t 03 

1 

583 

-8s« 

813 

—507 

587 

j -336 

393 

*51 

1 

AX * 

109 

*37 

—361 

116 

694 

-608 

• 54 

->6*9 


- 4*3 

457 

6 

xr 

XT 

39 * 

308 

-663 

4*1 

-341 

i *09 

— »®5 

109 

- 

- 5®3 

1 * 4 * 

-tsii 


! 

-891 

791 

tSs 

“^13 ^ 

1 —166 

I 

7*1 

—780 


! 

3 l >3 

— 4®9 

i 

^ AX * 

— 33 * 

— 6 

180 

•51 I 

1 — 61 ' 

1 -« 7 * * 

73 

—303 

438 

! -*78 

34* 

— 4tol 

lOOf* 

af 

- 5*9 

105 j 

-311- 

390 

1 — * 5 * 

! 153 

— 75 

431 

54 

j -*59 

753 

-*53 


XI 

—793 

-13* 

-141 

546 

j 

103 

j 400 

597 

- 45 

43 ® 

! ^ 

67 

3*6 

i 

' AX \ 

—*04 

i 

— 95 

1 

—114 

353 

* 3>4 ^ 

*78 

— 3 * 

554 

; - 4>5 

*70 

—103 

lUf 

XT 

*36 

—167 

*68 

3 * 

—168 

•54 

-*s« 

593 

1*4 

-319 

• 5 ® 

114 


; 

8*6 

• -537 

-639 

—499 

-794 

3*3 

i . 

6 

-351 

388 

1 3 * 

— 4 ■ 

43 ® 

! 

AX 

*47 

ns 

— *94 

—441 

-36 

' - 77 

45 ® 

—sot 

— S 3 

' 57 

— 4 

71 

laf 

AF 

—*380 

603 

641 

•4 

, —180 

*46 

- 33 ® 

>44 

— *»3 

- 15 * 

-336 

713 


AZ 

8*5 

14 * 

-ISS 

75 ® 

, -393 

— 44 

1 

-1101 

j -8S3 

i *»' 

• 

-354 

451 

— *74 

i 

i ^ 

396 

— 13 

— 97 

430 

•03 

j - 35 ® 

45 * 

-J.. 

j-m 

*88 

*63 

-1*8 

130 * ' 

! AF 

—174 

-783 

535 

-306 

—19* 

1 396 

—300 

! - 4*3 

t *51 

! 75 

—693 

615 


1 

4*5 

•94 

-387 

t 1613 

1 

i i ®54 

! 446 

4 

-ii6i 

lets 

I- 

•94 

. — 77 ® 

• 7 * 

—870 


i XX ■ 

81 

—*63 

45 

1 

< 

i 331 

' -*56 

*31 

-.75 

— 116 

5 «i 

— *®4 

— 3 » 

140 * : 

; XT : 

»9 

—341 

S «9 

1 — Sio 

1 -»47 

1 ’6* 

— 5 * 

1 —6*4 

40 

7 ®« 

—667 

iSo 


XI ‘ 

S6 

-.94 

—101 

j - 1*41 

1 - 3 P 9 ® 

: ’*♦ 

548 

*467 

439 

-489 

667 

-lOJD 


1 ^ 

—399 

1 

j -*03 

-613 

1 1*5 

1 — 7*4 

; —747 

394 

6*7 

677 

-639 

•34 

-478 

iXt 

XT . 

367 

! 419 

57 


j -180 

; 605 

- 6* 

-598 

6*6 

•771 

a*6 

35 * 



-83 

j-.S8 

74 * 

1 — 61 

! 1646 

: 768 

-•847 

>-1918 

—**31 

99 

* 4 *» 

«5 


*This table (not mentioned in the text) has been added to provide 
an overview for the next two tables (IX and X). It contains 
the differences of the computed and observed values of X, Y, Z, 
i.e. the values in VIII and those in ill, formed as (observed - 
computed). The corresponding differences for the intermediate 
points (» = 0°, 30°, 60°...330°) of the same parallel circles 
are found (according to provisional computation) in B, table 
XlVa, b, c, p. 64-66. 

















ORIGINAL PAGE li 
OF POOR QUALITY 


IX. Coefficients of the Series Representing 17, w, Vtk ft — a,»m. ) 

1 . 




m; n: 

0 

1 

2 

• 

4 

fi 

4 

0 

0 

-16428.t 

- 888.0 

854.4 

876.0 

- 58.0 

54 

1 


- 1167.0 
8465.T 

1902.0 
- S04.0 

-407.0 
— 92.T 

147.0 

61.T 

110.T 
— 894 

184 

884 

2 



808.0 

716.0 

587.4 

81.1 

189.1 
- 89.4 

106.4 
- 8.1 

- 74 
16.1 

8 




141.0 

885.0 

-100.0 
- 55j 

6.T 
— 04 

-854 

-104 

4 


i 



lOU 

15J 

1 1 
m •» 

54 

184 


17« Z7,+i7,(-633.»«*l + 60.o«nA) 4./7i(27jeM22»42.tMilS2) 
+ /7| (-202.J a»8i +154.t sin 8i)-M7« (89.»eot42-11.4«m41) 


IF. 


nr, ft: | 

0 

1 

2 

8 

4 

6 

6 

0 

0 

-16488.0 

- 888.0 

854.4 

8764 

- 584 

54 

1 


- 1277.0 
3481.0 

1861.0 
- 380.O 

-496.0 

-1324 

166.0 

1214 

1704 
- 614 

-884 

74 

8 



3414 

633.0 

585.0 
- 114 

186.0 
- 474 

1084 

0.0 

— 14 
884 

8 




172.7 

849.7 

-1044 
- 67.0 

8.0 
— 4.0 

-874 
— 5.7 

4 





404 
- 68.0 

— 214 

- 27.0 

184 

884 


W = (50- 67 J21-34j?J+38B*+9J?;+4.«7S;-0.tTiJ5) 


Vib 


m; ft: 1 0 

1 

2 

8 

4 

6 

6 

. 0~ 

-184284 

- 238.0 

8544 

8764 

- 584 

54 

1 

- 1222.0 

1831.0 

-4514 

1564 

1404 

— 94 

1 1 

84534 

- 317.0 

-1184 

914 

- 754 

19.0 

1 


275.0 

5864 

1874 

1074 

- 4.1 

8 1 

j 

6704 

44 

- 384 

— 84 

854 

1 



157.0 

-1084 

44 

-864 

3 1 



8374 

- 614 

- 8.0 

- 64 

1 




70.7 

- 88.0 

9.7 

* 1 




- 264 

- 864 

174 


84 



ORIGINAL PAGE rS 
OF POOR QUALITY 


IX. Coefficients of the Series Representing V, W,.7ik 0 >■ <JS6.10*cm.) 


2. 

U. 


m; n: 

1 ^ ' 

1 

2 

8 

4 

0 

0 

-18400.4 

- 287.0 

898.0 

271.1 



- 1701.1 

1254J 

-589.7 

156.0 

1 


8656.1 

— 888.7 

- 85.0 

48.0 




292.0 

488.0 

1853 

2 



C44.0 

7.7 

- 543 



■ ■ 


167.7 

- 89.7 

8 


' ' 


247.0 

- 663 






453 

4 





- 613 


n— 0«4-J7i(-397.4eMa - 8.t«mi)+l?i<-»7.T«*21-69.»itil2i) 


+/7»(-174.T«»8i+107A«i»8i) 


s. 


0 

f 

2 

0 

-186183 

- 898.1 


- 10193 
85173 

1180.1 
- 8783 



2003 

6403 


17— Z7,+I7i (-499.7 eM 2-28.0 Ml i) 




m; ft: 

0 

1 

2 

3 

0 

1 ^ 

-184003 

- 287.0 

893.0 

1 

! 

- 18123 
8452.0 

1250.O 
- 358.0 

-644.0 

193 

2 



2423 

6223 

478.0 
- 13 

8 



! 

167.7 

247.0 

4 I 







Ok 

1 

© 1 

1 

2 

271.1 0 

-186183 

- 3984 

.117.0 

— 1019.0 

10513 

503 

85173 

- 8963 

1953 

1 

2oa» 

- 50.0 


C40.0 


- 99.0 

— 60.T 


45.0 

61.0 

F — F,+2(58 -28.ro 


F- F,+2(49-642?J-345*+41.mB;+9.«iJ20 


V.h 


m; ft: : 

0 

1 

2 

8 

0 

0 

-18400.4 

- 287.0 

898.0 

1 


- 1507.0 
3554.0 

12523 
- 8483 

-591.0 
— 83 

2 



2673 

688.7 

4803 

83 

8 




167.7 

2473 

4 






V:h 


4 

0 

1 

2 

271.1 

0 

-186183 

- 898,1 

187.0 

463 

' 1 

- 10193 
8517.0 

11153 
— 8843 

1903 
- 523 



2003 

640.0 


- 94.4 

- 98.0 


45.0 
- 61.S 


* 


85 


































